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ABSTRACT 


A theory  of  three-dimensional  infinitesimal  isotropic  thermovisco- 
plasticity  is  proposed  and  investigated  in  several  homogeneous  deformations. 

The  theory  is  represented  by  both  a mechanical  constitutive  equation  and  a 
constitutive  assumption  for  the  heat  equation;  these  equations  are  separately 
postulated  but  are  coupled  through  their  common  linear  dependence  upon  the 
stress  rate  and  strain  rate  tensors  and  the  time  rate  of  temperature  change, 
and  through  their  common  nonlinear  dependence  upon  the  stress  and  strain 
tensors  and  the  absolute  temperature.  Throughout  this  theory  total  strain 
is  used  and  the  strain  tensor  is  not  decomposed  into  a sum  of  elastic  and 
inelastic  contributions.  The  concept  of  the  yield  surface  is  not  employed 
and  the  transition  from  initially  linear  thermoelastic  behavior  to  nonlinear 
inelastic  behavior  is  continuous. 

Through  qualitative  discussion  and  quantitative  examples,  this  proposed 
system  of  constitutive  equations  is  shown  to  represent:  initial  linear  thermo- 
elastic behavior  followed  by  inelastic  rate  sensitive  work  hardening; 
linear  thermoelastic  behavior  between  stress,  strain,  and  temperature  in 
pure  hydrostatic  loading;  initial  cooling  in  uniaxial  tension,  initial 
heating  in  uniaxial  compression,  and  initial  isothermal  behavior  in  torsion 
prior  to  inelastic  deformation;  subsequent  self-heating  in  any  monotonic 
loading  following  the  onset  of  inelastic  behavior;  axial  stress  (strain) 
build-up  due  to  temperature  change  uniformly  induced  by  monotonic  and  cyclic 
torsional  loading.  Also  represented  are  stress  (strain)  rate  sensitivity; 
temperature  rate  sensitivity;  creep  and  relaxation;  and  nonlinear  dependence 
of  the  spacing  between  stress-strain  curves  at  different  stress  (strain) 
rates  upon  the  value  of  the  stress  (strain)  rate  tensor  in  constant-rate 
loading.  Stress-strain  curves  corresponding  to  different  constant  strain 
(stress)  rate  loadings  ultimately  attain  the  same  slope.  Defined  and 
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physically  meaningful  limits  are  obtained  in  the  cases  of  very  slow  and  very 
fast  loading  rates  in  monotonic  radial  loading.  Associated  with  limitingly 
slow  loading  is  an  equilibrium  stress-strain  curve,  and  in  relaxation  the 
stress  relaxes  to  an  equilibrium  value  associated  with  this  curve. 

Illustrative  specializations  to  cases  of  uniaxial,  torsional,  pure 
hydrostatic,  and  radial  deformation  histories  are  considered  and  numerical 
results  are  obtained  by  postulating  specific  choices  of  material  functions, 
and  by  numerically  integrating  the  resulting  system  of  first-order  nonauto- 
nomous,  nonlinear  stiff  differential  equations.  Results  are  presented  in 
terms  of  graphs  to  show  the  influence  of  various  test  histories  and  material 
parameters  as  well  as  the  capability  of  the  model. 

In  cyclic  deformation  the  proposed  equations  of  the  model  are  modified 
according  to  previously  established  concepts  to  account  for  history  depend- 
ence in  the  sense  of  plasticity.  This  is  accomplished  by  updating  a material 
parameter  associated  with  the  equilibrium  stress-strain  curves 
in  such  a manner  that  the  material  response  predicted  by  the  model  remains 
continuous  and  smooth.  The  ability  to  represent  hardening  in  materials  is 
illustrated.  In  torsional  cycling  a significant  net  increase  in  temperature 
can  be  induced.  A possible  thermoinelastic  explanation  for  the  Ronay  effect 
is  given. 


INTRODUCTION 

Mechanical  deformation  of  materials  is  accompanied  by  deformation- 
induced  changes  in  temperature,  which  may  or  may  not  be  equilibrated  through 
heat  conduction.  Under  special  circumstances,  however,  such  as  rapid  mono- 
tonic and  cyclic  loading,  the  temperature  may  rise  considerably  and  may 
affect  material  deformation  behavior.  To  model  these  phenomena  a coupled 
mechanical  and  thermal  theory  of  material  behavior  is  necessary.  In  this 
paper  we  propose  a coupled  rate-dependent  isotropic  theory  in  the  context  of 
infinitesimal  deformation  based  on  total  strain  [1,2,3]. 

General  theories  of  thermomechanical  behavior  of  continua  have  been 
developed  in  [4-8] . A reduction  to  specific  classes  of  solids,  such  as  metals 
or  polymers  has  not  been  attempted.  Rate-independent  thermoplasticity 
theories  for  the  modelling  of  metal  deformation  behavior  were  proposed  in 
[9-12],  Ihermoviscoplasticity  was  considered  in  [13,14].  The  concept  of  a 
thermorheologically  simple  material  has  been  proposed  [15-19]  which  combines 
the  effects  of  temperature  and  time  into  one  parameter  of  process.  Other 
models  are  based  on  thermoviscoelasticity  [20-23]. 

Experimental  information  on  deformation-induced  temperature  changes 
is  not  abundant  [10-12, 24-29] . These  experiments  are  difficult  to  perform  and 
require  sophisticated  equipment  which  became  available  during  the  last  decade. 
Such  equipment  was  utilized  in  [30]. 

In  developing  the  coupled  theory  of  thermoviscoplasticity  based  on  total 
infinitesimal  strain  we  have  considered  the  experimental  results,  our  iso- 
thermal theory  [1-3]  and  the  first  law  of  thermodynamics.  This  theory  uses 
total  strain  only,  regards  rate  dependence  as  fundamental,  is  nonlinear  in 
the  stress  tensor,  the  strain  tensor  and  the  temperature,  but  is  linear  in 
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the  stress  and  strain  rate  tensors  as  well  as  in  the  temperature  rate.  The 
concept  of  a yield  surface  is  not  employed  and  the  transition  from  the  ini- 
tial thermoelastic  behavior  to  the  subsequent  thermoinelastic  behavior  is 
smooth. 

Ohe  mechanical  constitutive  equation  and  the  heat  conduction  are  de- 
veloped first.  Then  the  properties  of  the  proposed  equations  are  studied 
qualitatively  under  adiabatic  and  spatially  homogeneous  conditions.  Approp- 
riate limits  of  the  solutions  of  the  nonlinear  differential  equations  are 
obtained  for  limiting  loading  rates  (extremely  fast  and  extremely  slow)  and 
for  large  times.  Next  it  is  shown  that  instantaneous  large  changes  in 
loading  rate  result  in  a continuous  instantaneous  thermoelastic  response. 

The  equations  are  then  reduced  to  special,  spatially  homogeneous  deforma- 
tions, such  as  hydrostatic,  uniaxial,  biaxial,  and  torsional  states  of 
stress.  Numerical  experiments  using  hypothetical  but  realistic  material 
functions  and  constants  show  the  predictive  capability  of  the  model.  The 
case  of  torsional  cycling  is  of  special  interest  since  it  shows  that  tor- 
sional cycling  introduces  axial  displacements  or  stresses  depending  on  the 
boundary  conditions  due  to  deformation- induced  temperature  changes. 
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CONSTITUTIVE  EQUATIONS 
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The  Mechanical  Constitutive  Equation 

We  propose  this  constitutive  theory  within  the  context  of  infinitesimal 
deformation,  and  we  designate  e as  the  infinitesimal  strain  tensor  and  ct  as 
the  stress  tensor.  Hie  absolute  temperature  is  0 While  the  reference  tem- 
perature corresponding  to  a state  of  equilibrium  in  the  material  is  repre- 
sented by  0q.  A superposed  dot  denotes  differentiation  with  respect  to 
time,  and  the  usual  comma-notation  is  employed  to  represent  differentiation 
with  respect  to  position  coordinates.  Following  [1]  we  propose  the  iso- 
tropic thermomechanical  constitutive  equation 


m[CT,  e,  0]  i . . — k[c,  e,0]o.  . = ct.  . - G.  . [e,  0] 
Ti]  n li  n 


13  13  13 


with 


V€  6 . 

*ij  = (1  + v)  (1-  2v)  6ij  + <1  + v)  “ (1-  2v)  (e"  9o)6ij 


(1) 


(2) 
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In  the  above,  square  brackets  denote  "function  of",  a is  the  coeffi- 
cient of  thermal  expansion  Which  can  be  a function  of  temperature,  and  v is 
Poisson' s ratio  which  could  be  a function  of  the  stress  or  strain  tensors 
through  their  isotropic  invariants  and  a function  of  temperature.  These 
functions  must  be  obtained  from  suitable  experiments  or  may  be  approximated 
to  be  constant.  For  the  isotropic  theory  the  scalar  functions  m[  ] and  k[  ] 
depend  on  the  invariants  of  the  stress  and  strain  tensors.  They  are  re- 
stricted to  be  positive,  even,  bounded,  and  monotonic  functions. 

While  other  posibilities  exist  [1,2,3]  we  propose  to  make  m[  ] and 
k[  ] dependent  on  the  overstress  magnitude 

H]*5  (3) 


us-  z.v:--.w„rea?:a au.  ima,:. 
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and  we  propose  to  restrict  the  ratio* 


where  E is  the  modulus  of  elasticity.  These  specializations  are  particularly 
useful  for  modelling  metal  deformation  behavior  [1,2,3], 

Following  [2]  we  call  G the  equilibrium  stress-strain  curves  and  re- 
quire that  G be  odd  in  its  first  argument.  We  represent  the  function  G by 


writing 
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G.  .[« ,9]  = t.  • 
i]  ’ 10  9 
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where  the  symbol  9 denotes  a suitable,  isotropic  invariant  of  e and  tempera- 
ture taken  from  Table  1.  This  invariant  is  normalized  to  reduce  in  uniaxial 
deformation  to  the  uniaxial  strain,  and  as  an  example  of  a suitable  repre- 
sentation we  have  the  usual  "effective  strain" 

/l.  5 e . .e  . . . *5 

cp=cp  = ( (6) 

(1+V) 

where  e is  the  deviator ic  strain  tensor.  We  note  that  9 may  depend  upon 
both  j and  a (9  - 0Q)  6^.  The  function  g is  required  to  have  the  appearance 
of  a uniaxial  stress-strain  curve,  hence 

g [0,  0Q]  = 0 (7) 


^■aEiej  when  92-0.  (8) 

In  the  numerical  examples  we  use  a specific  representation  of  g[9,  9]  proposed 
in  [31]. 


Removing  this  restriction  is  discussed  in  the  isothermal  case  [3] . 


Equation  (1)  can  be  transformed  into  an  equivalent  integral  equation 


ds 


kf 


k[a[s] , 6 [s] , 9[s]  ] ! k l 


•V< 


» dT  (9) 


where  we  have  used  a[t=0]  =0,  e[t  = 0]  =0  and  0[t=O]  *0  , and  where  we 

<V  'V'  ~ O7 

employ  total  differentiation  with  respect  to  strain 

d = S d9 
dek*  3ekf  dekf  30  ' 

We  note  that  for  constant  material  properties 
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a*. . 


de 


v6.  -6w 

(6,1.6,  4 + 6,.^)  + 


d0 


kX  2 (1  + v)  ik  j£  ijJ  jk;  (1  + v)  (1  - 2V)  (l-2V)dek£ 


(10) 
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Thermomechanical  Coupling 

It  is  well  known  that  mechanical  working  causes  self-heating  of  ma- 
terials; in  some  cyclic  load  applications  inelastic  self-heating  can  result 
in  a temperature  high  enough  to  melt  the  material  [29] . A constitutive  theory 
should  be  able  to  represent  these  features  through  coupling  of  the  deformation 
with  the  equation  of  heat  conduction. 

The  deformation  of  metals  is  further  characterized  by  initial  thermo- 
elastic response  prior  to  inelastic  and  rate-sensitive  behavior.  During 
the  initial  elastic  response  the  deformation  produces  cooling  in  uniaxial 
tension  and  heating  in  uniaxial  compression,  While  pure  torsion  results  in 
no  temperature  change  (isothermal)  for  isotropic  materials.  Monotonic  deforma- 
tion into  the  inelastic  region  causes  only  self-heating.  Under  cyclic  load- 
ing the  temperature  fluctuates  around  a monotonically  increasing  mean 
temperature  [10-12,24,25,29,30], 

The  first  law  of  thermodynamics  for  infinitesimal  deformation  is 
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where  p is  the  constant  mass  density,  e is  the  internal  energy,  q is  the  out- 
ward heat  flux  vector,  q.  .is  the  divergence  of  the  heat  flux,  and  R is  an 
external  heat  supply  per  unit  mass. 

In  the  absence  of  external  heating  (R=0)  and  under  adiabatic  conditions 

(q.  . » 0)  the  change  in  internal  energy  is  entirely  due  to  the  mechanical 
1 

working,  i.e., 


p(e[t] 


-e[to])  = J a-.e..  dr. 


Thus  the  work  done  by  the  stresses  is  "stored"  in  the  material.  For  a load- 
ing starting  and  ending  at  zero  stress,  i.e.,  cr  ..  [t]  = [tQ]  =0,  two  forms 
of  storage  are  possible*;  change  in  temperature  and  change  in  the  microstruc- 
ture of  the  material  [28,32,33], 

Experimental  evidence  suggests  that  the  energy  of  microstructure  change 
is  a small  fraction  («  10%)  of  the  mechanical  working  [28,32,33]  and  for  the 
purposes  of  this  paper  we  assume  that  the  mechanical  working  is  completely 
transformed  into  heat  which  permits  us  to  calculate  the  temperature  increases 
due  to  self-heating  quite  well 

To  represent  this  behavior  we  postulate  that  the  internal  energy  rate 
is  determined  by 

pe-pC[61$  + & pgr  - 2f°L+<rtakk]  • <141 

where  C[6]  is  the  specific  heat  of  the  material  which  usually  is  a slowly 
varying  function  of  temperature. 


If  the  loading  would  start  or  end  at  nonzero  stress  then  mechanical  energy 
such  as  the  elastic  energy  could  be  stored  in  the  deformed  body. 

** 

Transformation  of  mechanical  working  into  microstructure  change  will  be 
considered  in  a future  paper. 


We  note  that  this  representation  of  the  energy  rate  is  an  exact  differ- 
ential and  this  permits  us  to  write 

e * efa(t),  6 (t) ;0  ] (15) 

o 

where  we  admit  the  parametric  dependence  of  energy  upon  the  initial  equi- 
librium temperature  and  where  we  set  e[O,0  ;0  ] = 0. 

7 o o 

Combining  the  first  law  of  thermodynamics  (12)  with  the  constitutive 
assumption  (14)  permits  us  to  write  our  equation  of  heat  conduction 

PCiejS-o  « - i o4Ja  - i tpK.  (16) 


It  is  seen  from  (14)  that  for  adiabatic  conditions  and  for  R=  0 the  mechanical 
work  in  a cycle  starting  and  ending  at  zero  stress  is  transformed  into  tem- 
perature change. 

Equations  (1)  or  (9)  and  (2),  (16)  represent  our  coupled  constitutive 
equations  of  thermoviscoplasticity.  These  follow  from  observations  of  ma- 
terial behavior  and  from  the  law  of  conservation  of  energy.  In  the  case  of 
nonadiabatic  deformation  we  must  additionally  postulate  a choice  for  the 
heat  flux  vector,  and  we  take  Fourier' s isotropic  heat  conduction  law  q^  = K0, i 
where  K is  the  thermal  conductivity. 


Hie  second  law  of  thermodynamics  may  be  written  as  [34] 


P0T]  i pR-q.  , + 


qi6.i 


Equation  (17)  does  not  determine  the  rate  of  entropy,  it  rather  specifies  a 
lower  bound  for  the  entropy  production  rate.  It  is  generally  understood  that 
this  entropy  inequality  must  not  be  violated.  We  select 

P9i|(t)  - PC (9) 6 . A (ilia  o a . i ) 


dt  \ 2E  ij  ij 
qi9  i 

” CTij®ij  + ~~t~  + A 


Equations  (1),  (2),  (16)  and  a suitable  law  for  the  heat  flux  vector 
constitute  a system  of  nonlinear  first-order  differential  equations  which 
are  linear  in  the  rates  of  the  stress  and  the  strain  tensors  and  the  temper- 
ature. For  their  solution  we  must  employ  the  equations  of  motion  and  we 
must  specify  the  external  heat  supply  and  the  boundary  conditions.  If  we 
consider  zero  heat  supply  and  spatially  homogeneous  conditions,  the  heat  flux 
vanishes  and  the  equations  of  equilibrium  are  satisfied.  We  then  must  either 
specify  the  stress  or  strain  tensor  as  a function  of  time.  The  temporal 
variations  of  the  strain  tensor  (if  the  stresses  are  prescribed)  or  the  stress 
tensor  (if  the  strains  are  prescribed)  and  the  temperature  are  then  obtained 
from  (1),  (2),  (16)  and  suitable  initial  conditions. 

For  q^  ^ = 0 and  R=  0 (16)  is  homogeneous  of  degree  one  in  the  rates  and 
can  be  expressed  as 


pc 


d0 


de 


ki 


d A + V „ 
ak i ' dekje  1 2E  CTijGij 


v 

2E 


(19) 


which  for  a = 0 reduces  to 

<v  /\/ 


de 


0f6  <jamm 
PC  dekf  ’ 


(20) 


For  constant  material  properties  (19)  can  be  rewritten  as 


1 + v ij  / V _ -A  ii 

°k l E ij  <3ekjg  \E  °mm  / de^ 

PC  + aCTkk 


Prom  (20)  we  see  that  the  model  predicts  initial  cooling  in  tension,  heating 
in  compression,  and  initial  isothermal  response  in  torsion  prior  to  the  on- 
set of  inelastic  response. 


Limiting  Behavior  for  Extremely  Slow  and  Fast  Loading 

We  impose  a constant  strain  rate  deformation  and  consider  the  limiting 
response  predicted  by  the  model  in  veiy  fast  and  very  slow  loading  at  con- 
stand  strain  rate.  Specifically  we  set 

®ij  * *ij  (22) 

with  e,  = e = f3t  and  Q, , = 1,  where  p and  Q are  constants. 

11  11  ^ 

Substitution  of  (22)  into  the  integral  equation  (9)  gives 

“ij  is. « - 1 ‘S. ' 0,  ♦ I { exp  - f 1 srr>  • ,23> 

6 

We  examine  the  response  under  infinitely  fast  (P -•<*>)  and  infinitely 
slow  (P-*0)  constant  strain  rate  deformation  by  performing  these  limits 


using  (23). 


We  obtain 


lim  a .=G  [e,0] 

P-0  J-j  ij  ~ 

CTi j " J Srfr  dT1  ^nn  d®  * 


Using  (4),  Eq. (25)  simplifies  to 


r d+n 

lira  ff.  . - E[0]  -r-±i  0 de 

ft-®  iD  t de 


and  to 


lim  o -Et 

0-»  «x>  ij  i] 


provided  E is  constant. 


Substitution  of  any  of  (24)  through  (27)  into  (21)  results  in  the 
appropriate  representation  of  the  heat  equation  corresponding  to  these 
limiting  cases.  From  (24)  we  identify  the  function  G as  the  equilibrium 
stress-strain  curves  of  the  material,  as  in  the  case  of  the  isothermal 
theory  [1] . We  see  that  (4)  does  not  guarantee  a linear  thermomechanical 
response  for  very  fast  deformation  unless  E is  independent  of  temperature . 

The  limiting  behavior  obtained  in  (24)  represents  a lower  bound  to  the 
mechanical  response  which  the  theory  predicts  at  any  finite  strain  rate, 
while  the  limits  (25)  - (27)  represent  am  upper  bound  to  this  response.  The 
same  limits  cam  be  obtained  for  infinitely  slow  and  infinitely  fast  constant 
stress  rate  loading.  These  limits  indicate  that  as  we  deform  the  material 
at  high  rates  the  response  tends  more  toward  elastic  behavior,  while 
deforming  the  material  at  low  rates  results  in  a response  closer  to  the 
equilibrium  response  manifested  in  G. 

Limiting  Behavior  at  Large  Times 


We  examine  the  behavior  predicted  by  this  theory  at  large  times  with 
constant  stress  rate  or  strain  rate  deformation.  Through  these  considerations 
we  obtain  limiting  bounds  upon  the  response  which  the  theory  predicts,  and 
while  these  limits  may  appear  to  imply  large  deformations,  they  are  actually 


rapidly  attained  asymptotic  limits.  Previous  investigations  [1-3]  and 
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present  numerical  experiments  indicate  that  these  limits  are  attained  by  the 
response  well  within  the  strain  magnitudes  appropriate  for  infinitesimal 
deformation  theory. 

We  define  the  difference  between  the  stress  and  the  equilibrium  re- 
sponse G to  be  the  overstress,  and  we  determine  the  limit  of  the  overstress 
as  time  goes  to  infinity  in  constant  strain  rate  deformation . We  use  the 
integral  equation  to  perform  this  limit  and  using  (4)  we  obtain,  following  [1,3] 

. df.  . dG.  .. 

Kj  - : : *'r>  “ (Ete)  - 1^)  \l  <28> 

where  the  braces  denote  the  time  limit  of  the  overstress 

fa.  . - G.  . } s lim  (a.  . - G.  . ) . (29) 

1 ID  ID  t-«  ^ a3 

This  limit  is  rapidly  attained  before  the  self-heating  temperature 
increase  is  sufficient  to  influence  the  mechanical  response  noticeably. 

Following  a similar  procedure  we  differentiate  (9)  in  time  obtaining 


the  limit  for  the  slopes 


t-»®  aejc£  t-*e 


dG, . [e,  0] 

lim  - A ~ = A 


ijk£  ‘ 


We  similarly  consider  the  limiting  overstress  in  constant  stress  rate 
deformation.  We  use  the  chain  rule  to  replace  the  strain  rate  term  in  the 
integrand  of  the  integral  equation  by  writing 


6.  . * -j— — a 
ki  da  mn 


We  proceed  by  taking  the  limit  of  the  overstress,  and  we  then  substitute  the 
"slope  limit"  (30)  into  this  limiting  overstress,  noting  that 
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Using  (32)  we  obtain  the  limiting  overstress  in  constant  stress  rate 


deformation 


, dt..  dG  . 

Kj-Gij5“kfr»0](E[ej  d^-di^)\^mn  * 


In  the  particular  case  of  our  representation  of  the  function  G proposed  in  [31] 


we  denote 


E [0]  = lim  flEPt91 
S cp-*«o  <P 


Then  we  have,  using  (34) 


A...  (-E  [61  3-=!  . 

iDkX  s deki. 


and  we  obtain  for  (33) 


taij"Gij}"kir’9l(Etei  _Es[01)  °ij/Es[01  * (36) 

Because  the  function  k[  ] depends  upon  the  special  invariant  defined 
in  (3),  Eqs. (28),  (33),  and  (36)  representing  the  limiting  overstress  are 
transcendental  equations,  indicating  that  the  limiting  overstress  is  finite 
and  depends  nonlinearly  upon  the  constant  strain  rate  or  constant  stress 
rate  tensors.  This  property  enables  the  modelling  of  arbitrarily  close 
spacing  of  stress-strain  curves  obtained  at  widely  different  stress  rates  or 
strain  rates,  as  shown  in  Fig.l  . If  the  function  k would  be  independent  of  T 
and  6 the  limiting  overstress  would  depend  linearly  on  the  applied  stress 
rates  or  strain  rates  [1,3]. 


Behavior  in  the  Neighborhood  of  the  Stress  and  Strain 
Origin  and  of  the  Initial  Temperature 


The  right-hand  side  of  Eq. (1)  is  initially  zero  in  an  undeformed  ma- 
terial, and  we  consequently  obtain  using  (4)  with  constant  material  properties 

d<7ij  d^ij 

aT1  " E aT1  * (37) 

deKjt  dekX 


Numerical  simulations  use  the  data  given  in  Table  2. 
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With  this  result  and  (20) 


2 

(1- 2v)pc- 3Ecr  0 


this  prediction  differs  from  classical  thermoelasticity  by  the  term  of 
2 

order  (Eor  0/pC) ; however  such  a term  is  much  smaller  than  one  (typically  less 

than  10  3)  and  its  effect  upon  the  prediction  is  negligible. 

2 

Consequently,  we  will  henceforth  neglect  the  term  in  our  inter- 

pretation of  the  heat  equation  and  we  regard  the  heat  equation  (38)  as  equiva 
lent  to  the  classical  thermoelastic  heat  equation. 


Instantaneous  Larqe  Chances  in  Strain  Rate  or  Stress  Rate 


Suppose  that  at  any  point  a,  e,  0 the  strain  rate  is  instantaneously 

changed  from  e..  to  eT.  where  e..  denotes  the  value  of  the  strain-rate  tensor 
13  13  lj 

prior  to  the  jump  change  in  strain  rate  and  eT..  denotes  the  value  of  the 
strain-rate  tensor  following  the  radial  jump  in  strain  rate.  We  are  inter- 

daij  dCTij  d0'  d0+ 

ested  in  the  change  in  the  slopes  -5 — to  — — - and  - to  at  the  instant 

dek  t dek  l dek£  dekJe 

the  strain  rate  is  changed.  From  (1)  and  the  chain  rule  we  obtain  using  (4) 
and  assuming  constant  material  properties 

. d^T.  da*.  , di|iT.  daT.v 

(e  - -T-^)  K,  “ (E  Ki  • (39) 

v dekJt  dek/  k*  v dek  1 dekr  k* 

To  represent  the  radial  jump  in  strain  rate  we  write  e . . ■ ae . . where 

ij  ij 

"a"  is  a constant,  and  obtain 


Although  it  is  lengthy  to  demonstrate  mathematically  the  approximations 
also  hold  when  the  material  properties  are  functions. 
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(40) 


Proceeding  similarly,  from  (21)  we  obtain 
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«ij(5  V-“9))  - \x(1  - I)}  • 


(41) 


2 *1 

If  | a|  » 1 (typical  changes  in  strain  rate  are  10  s and  higher)  we 


have  approximately 


4Si 


(42) 


and 


de 

de. 


oQE 


2 I 

kX  pc  (1  - 2v)  - 3er  0E 


(43) 


Comparing  (42)  with  (37)  and  (43)  with  (38),  respectively,  we  see  that 
large  radial  positive  or  negative  jumps  in  strain  rate  result  in  a thermo- 
elastic  slope  of  the  stress-strain  and  the  temperature-strain  diagram;  see 
Figs.l  and  2. 

To  consider  radial  jump  changes  in  stress  rate  we  invert  the  constitu- 
tive equation  (1)  with  (2)  and  using  (4)  in  the  case  of  constant  material 
properties  we  obtain 


eij  " Ek  [ ] (<1  + v)  (aij*Gij}  " v(akk“Gkk,6ij^  + E ((1+V)CT: 


ij 

kkvij/  “*6ij 


.)  + 006 


(44) 


• + 


We  set  o^  • bai_j  Where  "b"  is  a constant  and  obtain 


dgntn  = 1 fjn n . ( d9+  _ 1 d9~  \ / 1\ 

b a mnVdakjj  b d a^J  mnkA  to/ 


Equation  (16)  yields 


Then  for  | b|  » 1 

de+  .„+ 

rs  _ . d0 

dakj^  rskX  rs  dcr^ 

and 

d6+  ct9  . 
d5~  ="  PC  °ki  ‘ 


(45) 


(46) 


(47) 


(48) 


(49) 


The  results  in  (48)  and  (49)  are  analogous  to  (42)  and  (43),  respectively. 

Prom  this  analysis  we  see  that  large  changes  | b|  » 1 in  the  stress 
rate  produce  linear  thermoelastic  response.  We  note  that  this  important 
property  holds  irrespective  of  the  values  of  a,  e and  9. 


Reduction  of  the  Constitutive  Equations  to  Special 
Deformations  under  Adiabatic  Conditions 

In  examining  the  predictions  of  this  constitutive  theory  in  special 
deformations  we  asstime  for  simplicity  and  convenience  that  all  material 
properties  are  constant  unless  otherwise  specified.  The  predictions  in  the 
case  of  nonconstant  material  parameters  follows  directly  from  the  theory 
but  is  cumbersome  in  its  representation. 


- 


Pure  Hydrostatic  Deformation 

In  this  case  we  have  <7^-06^  and  e „ - e 6^.  Depending  on  the  choice 
of  the  invariant  cp  in  G different  behaviors  can  be  modelled.  If  a devia- 
toric  cp  from  Table  1 is  chosen  then  from  (1)  for  constant  material 


properties 


and  from  (21) 


• E*ij 


aOE . 

2 

(1- 2V)pC- 3E«  9 


and  we  see  that  the  behavior  is  linear  thermoelastic.  If  we  chose  a non- 


deviator ic  invariant  then  nonlinear,  rate-dependent  behavior  will  be  obtained. 
No  information  regarding  actual  material  behavior  under  nonisothermal 


hydrostatic  conditions  is  available  and  in  keeping  with  tradition  we  choose 
a deviatoric  invariant  for  metals.  For  other  materials  a nondeviatoric  cp 
appears  to  be  suitable. 

Phiaxial  Deformation 

In  this  case  - x(t)  - e - or(0  - 0q)  and  all  other  f„-0»  a^“ff(t) 
and  all  other  a^-0.  Equations  (1),  (2),  (6)  and  (16)  reduce  to 


a . 22  (|ta6) 

d0  m dc  'E  - 

d«  pc  + ota 


or  alternatively  we  write  (53)  as 


(pC-Bcr  0)0-  (a  + Ecr0)  (a-g[x,0J) 
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The  uniaxial  stress  response  and  temperature  response  under  constant  strain 
rate  are  shown  in  Figs. 3 and  4 (see  Table  2 for  data). 


Torsion  of  a Thin-Walled  Tube 

In  pure  torsion  of  a tube  one  usually  specifies  the  shear  stresses  or 

shear  strains  and  in  this  case  we  have  Y ■-  and  ct„  = <r  = t. 

12  21  1+V  12  21 

Inelastic  torsional  deformation  results  in  a temperature  increase  [10-12,24-28], 
and  this  temperature  change  induces  an  axial  "thermal  strain".  Consequently, 
in  torsion  we  must  also  examine  the  axial  component  of  the  constitutive 
equation  (1)  under  specified  boundary  conditions.  From  (1)  - (4)  and  (16) 
we  have 


& * x 

(i  +v)  1 T - krr,e] 

(T-S12), 

(56) 

/ C7\ 

E(e  ore)  a kfr)ej 

w / / 

Q ~ Q 

2a[t]>  T(T-S12>-E”9'  + EklV  + 

(58) 

In  the  above  (5)  reduces  to 


G - (e  - cr(6  - 0 ) ?[tP»ei 
11  o cp 


„ _Y_  SLLSEiil 

^12  1+v  9 


We  have  two  choices  for  the  axial  boundary  conditions  s zero  stress  or 
zero  strain.  If  we  require  zero  axial  stress  then  o and  cr  vanish  for  all  time 


Specification  of  nonzero  axial  stress  or  strain  requires  prior  axial  loading 
which  we  do  not  consider  here. 
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and  this  requires  ■ 0 and  G^  = 0.  Consequently  under  this  boundary  condi- 
tion the  axial  deformation  has  the  simple  solution 

e = or©  . (61) 

Equation  (57)  is  then  trivially  satisfied  and  Eq. (58)  reduces  to  the  torsional 
heat  equation  and  contains  only  shear  contributions 

■£#ffT<T-=12)  • (62) 

Thus,  when  assuming  this  boundary  condition  the  theory  predicts  a uniform 
elongation  equal  to  <*(0-0^). 


If  we  alternatively  assume  the  zero  strain  boundary  condition  (e  = 0, 
e = 0)  (57)  reduces  to 


a + 


a 

k[  ] 


9 [cp.  61  > 
k[  ] / 


(63) 


while  in  (58)  the  e-term  vanishes.  Equation  (58)  predicts  a temperature  in- 
crease due  to  the  inelastic  torsional  deformation  represented  by  (56).  This 
temperature  increase  results  in  an  axial  compressive  stress  build-up  according 
to  (63)  and  this  axial  deformation  can  further  contribute  to  the  self-heating 
of  the  material  as  is  clear  by  comparing  (58)  with  (62). 

In  the  first  case  the  axial  unit  elongation  is  proportional  to  the 
temperature  change  [see  (61)],  Which  is  generated  through  the  torsional  de- 
formation through  (62).  However,  in  the  second  case  there  is  zero  axial  unit 
elongation  but  there  is  a rate-dependent,  inelastic,  and  gradual  compressive 
stress  build-up  predicted  by  (63). 

Figures  1,  2,  3 and  4 show  the  stress-strain  and  temperature-strain 
response  under  pure  torsional  deformation  of  a thin-walled  tube  with  boundary 
conditions  of  zero  axial  stress  [Eqs.(56),  (61)  and  (62)  apply]. 


H 
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Strain-Controlled  Torsional  Cycling 

It  is  a fundamental  assertion  of  our  theory  that  viscoplastic  behavior 
can  be  represented  by  discontinuously  accounting  for  the  deformation- induced 
microstructure  change  of  the  material  [38,39].  In  our  theory  this  accounting 
is  accomplished  at  physically  defined  points.  The  values  of  stress,  strain, 
and  temperature  at  these  points  are  stored  in  material  functions  and  simul- 
taneously appropriate  constant  material  parameters  are  updated.  Storage  with- 
out updating  occurs  in  a cyclic  steady  state  condition.  Without  this  storage 
and  updating  our  equations  are  of  a nonlinear  thermoviscoelastic  nature. 

Through  this  procedure  we  represent  thermoviscoplasticity  with  piecewise 
nonlinear  thermoviscoelasticity. 

Although  storing  and  updating  is  discontinuous  the  response  predicted 
by  our  theory  remains  continuous.  We  specify  that  storing  and  updating  occur 
when  any  component  of  the  overstress  vanishes  and  subsequently  changes  sign. 

The  rules  for  updating  must  be  obtained  from  cyclic  deformation  experiments. 
Presently  the  rules  for  storing  and  updating  are  only  developed  for  cyclic 
proportional  loading. 

To  illustrate  the  procedure  we  simulate  in  Fig. 5 cyclic  hardening  in 
completely  reversed  strain  controlled  torsional  cycling  of  an  axially  con- 
strained thin-walled  tube  (e=*0). 

Equations  (56),  (58)  (with  e ■ 0)  and  (63)  apply  and  the  equilibrium 
functions  are  given  by  (59)  and  (60).  At  the  first  occurrence  of  vanishing 
torsional  overstress  the  values  of  the  equilibrium  functions,  strain,  and 
temperature  are  stored.  Simultaneously  the  physically  defined  parameter  Xf 
used  in  the  representation  of  g (see  [31])  is  updated*  to  simulate  hardening**. 

♦ 

Additional  possibilities  exist  through  the  updating  of  k[  ]. 

** 

If  necessary  more  than  one  parameter  can  be  changed  (see  [40] ) . 
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At  the  point  of  vanishing  torsional  overstress  we  denote  the  axial  equilibrium 
function  value,  the  shear  equilibrium  function  value,  the  shear  strain  and 
the  temperature  as  G.^,  y ^ and  0^.  Departing  the  point  of  van- 

ishing torsional  overstress  we  have 


V - Y. 


(1) 


12 


1 + V 


g(1)  [cp(1).ei 


9 


(1) 


+ G 


(1) 

12 


and 


(64) 


G^-cs-e11’)  .<;»>  (65 

9 

where  cp^  is  determined  using  y*Y^  and  g^  represents  the 

updated  stress-strain  curve. 

At  the  subsequent  occurrences  of  vanishing  torsional  overstress  the 
storage-updating  procedure  is  repeated. 

The  cyclic  deformation  shown  in  Fig. 5 results  in  the  temperature  in- 
crease depicted  in  Fig. 6 which  in  turn  leads  to  the  compressive  stress 
build-up  in  Fig. 7. 


DISCUSSION  AND  NUMERICAL  EXPERIMENTS 

The  preceding  constitutive  theory  represented  by  (1),  (2),  (16)  with 
(3)  reproduces  in  a unified  way  many  of  the  qualitative  features  of  non- 
isothermal,  inelastic  metal  deformation  generally  attributed  to  thermovisco- 
plastic behavior.  Ihe  capabilities  of  the  theory  include: 

• Initial  thermoelastic  response  at  all  loading  rates  prior  to 
the  onset  of  inelastic  behavior. 

• Subsequent  deformation-induced  temperature  increase  during  the 
nonlinear  (inelastic)  mechanical  response. 

• Linear  thermoelastic  behavior  at  all  values  of  stress  and 
strain  in  the  pure  hydrostatic  deformation  field. 


Stress  rate,  strain  rate  and  temperature  rate  sensitivity 
of  the  mechanical  response. 

Creep  and  relaxation  are  treated  in  a unified  way. 

Stress  rate  and  strain  rate  sensitivity  of  the  thermal  response. 

Initial  thermoelastic  "slope"  of  both  mechanical  response  and 
thermal  response  immediately  upon  the  imposition  of  large 
instantaneous  increases  in  the  strain  rate  or  stress  rate. 

Defined  mechanical  and  thermal  behavior  under  limitingly  slow 
and  limitingly  fast  rates  of  deformation. 

Nonlinear  dependence  of  steady-state  mechanical  response  upon 
the  applied  constant  stress  rate  or  constant  strain  rate 
tensors  (nonlinear  spacing  between  stress-strain  curves). 

Absence  of  strain  rate  history  effect  on  stress;  presence  of 
strain  rate  history  effect  on  temperature. 

Stress-strain  curves  obtained  at  different  constant  strain  rates 
will  ultimately  have  the  same  " slope" . 

Temperature  dependent  "softening"  of  both  initial  elastic  and 
subsequent  inelastic  response. 

Significant  net  increase  of  temperature  under  adiabatic 
cycling  in  torsion. 

Cycling  hardening.  Permanent  set  upon  unloading  to  zero  stress. 

In  loading  from  and  unloading  to  zero  stress,  energy  of  mechanical 
working  is  completely  converted  to  temperature  change. 

Temperature  induced  coupling  of  torsional  and  axial  deformation 
as  a possible  source  for  the  Ronay  effect.  Torsional  cycling 
is  demonstrated  to  result  in  temperature  increase  which  iso- 
tropically induces  either  an  axial  strain  field  or  an  axial 
stress  field. 


In  the  above  the  terms  "slope"  and  stress-strain  refer  to  particular 
components  of  tensors  specialized  to  one-dimensional  homogeneous  deformations 
such  as  the  tensile  or  torsional  tests. 

Ihe  above  capabilities  are  substantiated  by  the  theory  and  are  illus- 
trated by  numerical  simulations.  In  the  numerical  examples  we  select  material 
functions  and  integrate  the  resulting  coupled  nonlinear  stiff  differential 
equations  using  a computer  program  based  on  an  algorithm  developed  in  [41] . 

The  material  functions  do  not  represent  a specific  real  material.  They 
are  listed  in  Table  2 and  are  chosen  to  approximate  the  behavior  of  a low 
strength  steel.  Unless  otherwise  mentioned  these  functions  are  used  in  the 
computation  of  the  responses  shown  in  the  following  figures. 

In  Figure  1 we  simulate  pure  torsion  at  constant  shear  strain  rates. 

The  capability  of  nonlinear  spacing  in  the  model  is  clear,  since  increasing 
the  strain  rate  two  and  four  orders  of  magnitude  in  this  simulation  produces 
a small  increase  of  steady-state  flow  stress.  This  behavior  is  predicted  under 
the  limits  (28)  and  (33).  Further,  all  stress-strain  curves  become  parallel 
to  each  other  and  to  the  equilibrium  function  G12,  as  predicted  by  the 
limits  (30)  and  (35).  At  1%  shear  strain  an  instantaneous  jump  in  shear 
strain  rate  is  imposed.  Ihe  stress  response  immediately  departs  the  point 
of  the  jump  in  rate  with  a thermoelastic  slope,  as  predicted  by  Eqs.  (42),  (43). 

We  also  see  in  Figure  1 and  from  the  theory  that  there  is  no  strain  rate 

history  effect  on  stress  since  the  response  corresponding  to  the  jump  increase 

rapidly  merges  from  below  with  the  response  corresponding  to  uniform  applica- 

tion of  the  high  strain  rate. 

In  Figure  2 we  have  the  temperature  history  corresponding  to  the  shear 
deformations  of  Figure  1.  We  see  in  Figure  1 that  deformation  at  increasing 
constant  strain  rates  increases  the  region  of  apparent  initial  linear 


mechanical  response.  During  the  initial  linear  mechanical  response  the 
thermal  response  predicted  by  (16)  or  (21)  is  thermoelastic  and  consequently 
Figure  2 shows  isothermal  response.  Subsequently  the  inelastic  torsional 
response  results  in  self-heating,  and  deformation  at  increasing  strain  rates 
results  in  increased  heating.  In  Figure  2 we  see  the  temperature  history 
corresponding  to  the  jump  in  strain  rate.  While  the  stress  response  has  no 
strain  rate  history  effect,  the  thermal  response  has  a clear  strain  rate 
history  effect.  Ibis  strain  rate  history  effect  on  temperature  response 
allows  distinction  between  the  continuous  and  jump  histories. 

Figure  3 indicates  the  predictions  of  this  constitutive  theory  in 
tension,  compression,  torsion,  and  proportional  loading.  In  Figure  4 we 
have  the  corresponding  temperature  responses  predicted  by  (16).  While  the 
temperature  and  temperature  rate  are  higher  in  compression  than  in  tension 
(Fig. 4)  the  mechanical  response  in  this  case  is  negligibly  affected  and  the 
two  stress  responses  coincide  on  the  graph  in  Figure  3.  We  see  that  the 
initial  response  is  thermoelastic  but  when  the  mechanical  behavior  becomes 
inelastic,  self-heating  rapidly  occurs.  In  the  case  of  proportional  loading 
the  mechanical  response  is  decreased  while  the  inelastic  thermal  response  is 
increased  relative  to  the  pure  shear  and  pure  tensile  response. 

In  Figure  5 we  simulate  torsional  cycling  of  a hypothetical  hardening 
material  using  shear  strain  control;  also  shown  is  the  torsional  equilibrium 
function.  When  the  torsional  overstress  vanishes  we  execute  the  updating 
procedure  described  previously.  To  simulate  hardening  the  parameter  X^ 
continues  to  change,  see  Table  2.  In  a cyclic  steady  state  (closed  loop)  the 
value  of  X^  is  unchanged  although  storage  continues.  We  note  that  the  response 
remains  smooth  even  at  the  update  points.  If  cycling  would  be  stopped  at 
zero  stress,  an  aftereffect  would  occur.  The  strain  magnitude  would  decrease 
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until  reaching  a nontrivial  equilibrium  value  determined  by  the  intersec- 
tion of  the  corresponding  equilibrium  curve  with  the  strain  axis.  A perm- 
anent set  is  predicted  by  our  model. 

In  Figure  6 we  see  the  adiabatic  temperature  response  corresponding 
to  the  torsional  cycling  in  Figure  5.  We  note  the  presence  of  both  heating 
and  isothermal  phases  during  each  torsional  cycle  similar  to  the  experimental 
curves  reported  in  [10-12,  24,25]. 

Tfais  temperature  response  gives  rise  to  a comprehensive  stress  build- 
up since  the  axial  elongation  was  required  to  be  zero.  This  axial  stress 
build-up  is  shown  in  Figure  7.  While  the  magnitude  of  the  axial  stress  is 
well  within  the  linear  region  of  Figure  3,  the  behavior  depicted  in 
Figure  7 is  inelastic.  In  the  case  of  zero  axial  stress  an  axial  elongation 
occurs  according  to  Eq.  (61).  This  thermally- induced  behavior  could  be  an 
explanation  of  the  Ronay  Effect  [35-37] . 

In  Figure  8 we  have  pure  tensile  response  at  one  constant  tensile 
-3  -1 

strain  rate  (10  s ) but  at  different  initial  temperatures.  The  effect 
of  "softening"  the  material  arises  through  decrease  of  the  elastic 
modulus  E[0]  with  increasing  temperature  depicted  in  the  insert..  In 
Figure  9 we  have  the  corresponding  temperature  response.  Increasing 
initial  temperature  results  in  increasing  intervals  of  elastic  cooling 
prior  to  inelastic  temperature  increase  as  predicted  by  (16),  (53)  or  (54). 

Examination  of  Figures  1-4  and  8,  9 shows  that  the  temperature 
increase  due  to  self-heating  is  usually  small  and  of  negligible  influence 
on  mechanical  behavior.  However,  extensive  adiabatic  cycling  may  increase 
the  temperature  considerably  with  significant  influence  on  mechanical 
behavior. 
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me  equations  proposed  herein  are  linear  in  the  temperature  rate.  In 
case  of  external  heating  the  rate  of  heating  will  have  an  influence  on  the 
mechanical  deformation.  Such  influences  of  temperature  rate  on  mechanical 
deformation  are  reported  in  [42]  for  the  case  of  external  heating  and 
mechanical  deformation  of  a uniaxial  bar.  Equation  (52)  shows  that  both 
the  rate  of  straining  and  the  temperature  rate  have  an  effect  on  the  stress 
It  predicts  an  increase  in  stress  with  increasing  temperature  rate  as  ob- 
served in  [42] . In  this  case  the  heat  source  term  R must  be  added  to 
Eq. (53). 
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MATERIAL  FUNCTIONS  AND  PROPERTIES  USED  IN  THE  NUMERICAL  EXPERIMENTS 


9 f<P,  6Jm  ES<P  + 2R  (tanh  {RX  _ 3 , ) lo9e(cosh  (V) ) 


In  torsional  cycling,  Figures  5-7,  the  parameter  Xf  is  updated  at 
successive  occurrences  of  vanishing  overstress.  Successive  values  of 
X„  are:  .0030,  .0060,  .0065,  .0069,  .00725,  .00745,  .00760,  .00765. 


In  cycling  this  parameter  value  is  increased  in  the  update- storage 
procedure. 
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Figure  8 

Figure  9 


Shear  Response  under  Constant  Shear  Strain  Rates. 

The  equilibrium  stress-strain  curve  G12  and  the  elastic 

response  are  also  shown  together  with  a response  to  an 
instantaneous  change  in  strain  rate.  The  merging  of 
curves  4 and  5,6  indicates  absence  of  strain  rate  history 
effect  on  stress. 

Adiabatic  Temperature  Response  for  the  Constant  Strain 
Rate  Tests  2 - 6 in  Fig. 1.  Curves  4 and  5,6  never  merge 
indicating  presence  of  strain  rate  history  effect  on 
temperature. 

Stress-Strain  Response  (Absolute  Values)  for  Tension, 
Compression,  Shear,  and  Combined  Proportional  Loading 
with  e a Y. 

Adiabatic  Temperature  Change  in  °K  vs.  Absolute  Strain 
for  the  Tests  Depicted  in  Fig. 3. 

Torsional  Strain  Cycling  of  a Thin-Walled  Tube  Under 
Axial  Constraint  (e  = 0)  at  y=*012  sin  2rrt*  Cyclic 
hardening  is  simulated  and  corresponding  equilibrium 
curves  G^  316  shown. 

Adiabatic  Temperature  Change  Corresponding  to  Torsional 
Cycling  Shown  in  Fig. 5. 

Compressive  Axial  Stress  Increase  in  the  Thin-Walled 
Tube  due  to  Torsional  Cycling.  The  equilibrium  curve  G^ 
is  also  shown. 

The  Influence  of  Initial  Temperature  and  Temperature 
Dependent  Elastic  Modulus  on  the  Stress-Strain  Diagram 
in  Tension,  e ■ 10"3  s“^.  The  temperature  dependence 
of  the  elastic  modulus  is  shown  in  the  insert. 

The  Temperature  Change  Corresponding  to  the  Deformations 
Shown  in  Fig. 8. 


Stress  response 
Corresponding  G. 
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absolute  temperature,  throughout  this  theory  total  strain  is  used  and  the 
strain  tensor  is  not  decomposed  intb  a sum  of  elastic  and  inelastic  contribu- 
tions. the  concept  of  the  yield  surface  is  not  employed  and  the  transition 
from  initially  linear  thermoelastic  behavior  to  nonlinear  inelastic  behavior 
is  continuous.  L_s—*** 

through  qualitative  discussion  and  quantitative  examples,  this  proposed 
system  of  constitutive  equations  is  shown  to  represent:  initial  linear  thermo- 
elastic behavior  followed  by  inelastic  rate  sensitive  work-hardening;  linear 
thermoelastic  behaviok  between  stress,  strain,  and  temperature  in  pure  hydro- 
static loading;  initials  cooling  in  uniaxial  tension,  initial  heating  in  uni- 
axial compression,  and  initial  isothermal  behavior  in  torsion  prior  to  inelastii 
deformation;  subsequent  stelf-heating  in  any  monotonic  loading  following  the  on- 
set of  inelastic  behavior; -axial  stress  (strain)  build-up  due  to  temperature 
change  uniformly  induced  by ^monotonic  and  cyclic  torsional  loading.  Also  rep- 
resented are  stress  (strain) \pate  sensitivity;  temperature  rate  sensitivity; 
creep  and  relaxation;  and  nonlinear  dependence  of  the  spacing  between  stress- 
strain  curves  at  different  stress  (strain)  rates  upon  the  value  of  the  stress 
(strain)  rate  tensor  in  constant-rate  loading.  Stress-strain  curves  correspond- 
ing to  different  constant  strainX (stress)  rate  loadings  ultimately  attain  the 
same  slope.  Defined  and  physically  meaningful  limits  are  obtained  in  the  cases 
of  very  slow  and  very  fast  loading -rates  in  monotonic  radial  loading.  Associ- 
ated with  limitingly  slow  loading  is  an  equilibrium  stress-strain  curve,  and  in 
relaxation  the  stress  relaxes  to  an  equilibrium  value  associated  with  this 
curve. 


Illustrative  specializations  to  closes  of  uniaxial,  torsional,  pure  hydro- 
static, and  radial  deformation  histories  are  considered  and  numerical  results 
are  obtained  by  postulating  specific  choices  of  material  functions,  and  by 
numerically  integrating  the  resulting  system  of  first-order  nonautonomous,  non- 
linear stiff  differential  equations.  Results  are  presented  in  terms  of  graphs 
to  show  the  influence  of  various  test  histories  and  material  parameters  as  well 
as  the  capability  of  the  model.  \ 

In  cyclic  deformation  the  proposed  equations  of  the  model  are  modified 
according  to  previously  established  concepts'^ to  account  for  history  dependence 
in  the  sense  of  plasticity.  This  is  accomplished  by  updating  a material  param- 
eter associated  with  the  equilibrium  stress-strain  curves  in  such  a manner  that 
the  material  response  predicted  by  the  model  remains  continuous  and  smooth. 

The  ability  to  represent  hardening  in  materials  is  illustrated.  In  torsional 
cycling  a significant  net  increase  in  temperature  can  be  induced.  A possible 
thermoinelastic  explanation  for  the  Ronay  effect  is  given. 
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